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ABSTRACT 

In the Standard Model (SM), dilepton production in hadron-hadron col- 
lisions proceeds through the conventional Drell-Yan mechanism qq — > 
with the exchange of a gauge boson. Some extensions of the SM contain 
a quark-lepton contact interaction via a qlip Yukawa coupling, where is 
a scalar. Theories with scalar leptoquarks and i?-parity violating SUSY 
models are the most important examples of such extensions. These Yukawa 
couplings induce a different dynamical configuration compared to the SM 
(t-channel vs. s-channel) in the qq — > l^l^ process and thus offer the possi- 
bility of being identifiable upon imposition of suitable kinematic cuts. We 
discuss these effects in the context of the dilepton production in the CDF 
experiment, and explore consequences in the forthcoming Large Hadron 
Colhder (LHC). 
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In the quest of understanding physics at energy scales beyond that of the standard 
model (SM), many ideas have been discussed rather extensively in the literature. Two 
of the most popular, and in some sense intertwined, theories are those of grand uni- 
fication and supersymmetry. In both theories, there arise new particles that connect 
the quark and lepton sectors. Such a particle can either be a scalar or a vector. In 
this letter, we shall restrict ourselves to scalars. As for weak SU (2) properties, it can 
transform either as a singlet, a doublet or a triplet. 

In grand unified theories [|l] , where quarks and leptons are part of the same multi- 
plet, there naturally arise scalar leptoquarks leading to a contact interaction between a 
quark and a lepton. In the Minimal Supersymmetric Standard Model (MSSM) too, 
the squarks may mediate interactions between quarks and leptons unless some extra 
symmetry (i?-parity) is imposed on the theory. As there is no deep motivation for such 
a symmetry to actually exist, it is interesting to examine the possible consequences of 
such i?-parity violating (-^p) interaction. Within the context of the MSSM, the part 
of the Lagrangian violating i?-parity may be written in terms of the chiral superfields 
as 

^fi = KjkLiLjEl + X[ji^LiQjDl + X'-ji^UiDjDl, (1) 

where Li and Qi are the SU{2) doublet lepton and quark fields and E^,Uf,D^ are 
the singlet superfields. Concentrating on the X[j^ piece above, it is obvious that the 
phenomenology due to such terms is very similar to that of certain scalar leptoquarks. 

A generic scalar coupling between a lepton (/) and a quark (g) can be parametrized 

as 

C = glih^PL + hRPR)qct) + gHh'^PL + h'^PR)qct)' + h.c, (2) 

where g is the weak gauge coupling constant, is the charge conjugate field and 
Pl,r = (1 =F 75)/2 (0 and 0' are independent fields). Since diquark couplings violate 
baryon number, non-observation of proton decay forbids the simultaneous presence of 
diquark couplings as welLQ Furthermore, to suppress flavour- changing neutral current 
processes |Q, each such scalar is assumed to couple primarily to only one family of 
quarks (leptons) . Finally, to survive stringent bounds coming from helicity suppressed 
n ^ eu and similar processes Q], the scalar coupling is assumed to be chiral {i.e. it 
is not allowed to simultaneously couple to the left- and right-handed quarks)^. There 
is a subtlety in the case of the MSSM as the squarks fi and fn will always mix. One 
may argue that such effects are proportional to mj /msusY and hence would be small 
in most cases including the one that we are about to discuss. Still, the existence of 

-'^This implies that in the case of I^p theories A" cannot simultaneously exist with either of A or A'. 
However, A and A" are not relevant to the rest of our discussion. 

^In principle, couplings of scalars with the left-handed quark doublets cannot be made fully flavour 
diagonal in the quark sector, since the presence of CKM mixings prevents its coupling to be diagonal 
in the up- and down-sectors simultaneously. 



1 



this mixing does affect some rare processes, and, as a result, the bounds on I/lp 
couphngs are somewhat stronger than those for a chirally couphng leptoquark. With 
the above assumptions then, there are only five possible scalar couplings involving a 
charged lepton and these are listed in Table |l|. We have grouped scalars with different 
SU{2)l ® U{X)y properties together as their gauge interactions are irrelevant for the 
processes that we intend to study. 

Leptoquark couplings as well as I/ip couplings are constrained by direct searches 
at various colliders 0. The Fermilab CDF and DO experiments rule out a 
'first generation scalar' upto ~ 100 GeV, almost irrespective of the size of the Yukawa 
couplings. On the other hand, the HERA collaboration sets a limit of ~ 180 GeV 
1^ on the mass of a scalar coupling with electroweak strength to a first generation 
quark and electron. Similar bounds also hold for the case. Due to the presence of 



other decay channels for the squarks, the exact constraints are somewhat different [jTO 



For example, the Tevatron dilepton [|TT| data have been used to set lower limits on 
squark/gluino masses for almost any value of the $,p coupling. For scalars coupling to 
the third generation, the bounds are understandably weaker. The strongest constraints 
to date can be inferred from the loop effects on the leptonic partial widths of the Z 
12| . Stronger constraints on such scalars can be obtained in future colliders by looking 
either at r-number violating processes at LEP200 or at tt production |T^ at the 
Next Linear Collider. 

It is obvious that such experiments (other than those at HERA) have very little 
to say about the Yukawa couplings per. se. In this letter, we point out an experiment 
that would do so. To do this, we investigate the possibility of identifying the effects of 
a virtual scalar exchange in dilepton production in hadronic reactions. The signature 
we focus on is a lepton pair without any missing transverse momentum. The lowest 
order SM process leading to such a final state is the Drell-Yan mechanism i.e. qq — > 
(7*, Z*) l^l^ . The presence of such a scalar would introduce an additional t-channel 



diagram, an obvious consequence being a modification of the total cross section |T5 
A more sensitive probe, though, could be a comparison of the differential distributions. 
Owing to the extra diagram being a t-channel one in contrast to the s-channel SM 
contribution, it is conceivable that the effect of the former would be more pronounced 
for certain phase space configurations. 

A particular distribution viz. the rapidity dependence suggests itself. In the sub- 
process centre-of-mass (cm.) frame, one expects the t-channel contribution to be more 
central than the s-channel one. Another differential variable of interest could be the 
invariant mass of the lepton pair. Indeed, the CDF collaboration has studied the latter 



distribution for both dielectron and dimuon pairs [|T^. In the first part of this letter, 
we show how this experiment can discriminate between the SM and a theory with an 
extra scalar interaction. We compare the theoretical predictions with the existing data 
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and point out the improvements necessary to make our conclusions more quantitative. 
In the second part, we examine the same effect for the LHC, and speculate on the 
sensitivity that can be achieved. 

In principle, the contribution of a scalar of the type 0' (see eq. (0)) can be distin- 
guished from that of type only if (i) the colliding beams are asymmetric and (ii) the 
individual lepton charges are distinguished. We shall not dwell upon this possibility 
here. We make the simplifying assumption that at best one of the scalars in Table |l] is 
present. As an example, we present here the differential cross section for the process 
qq l^l^ for the case = 0, ^ 0. This corresponds to scalars of Types III, IV 
or V in Table 1. (The expressions for the other types of scalars can be obtained by 
simple modifications of the couplings.) Defining the left- and right-handed couplings 
of the Z to a fermion (/) by 



(tl — Cf sin^ 9w)/ sin 9w cos9w, 



—Cf cot 6w 



(3) 



where Cf and 4 are the corresponding charge and isospin respectively and 6*14^ is the 
weak mixing angle, we have, for a scalar of mass m<^: 



da 
dt 



dcT' 



SM 



dt 



3M6 



(M2 - m|)2 + r|m| ^ 4 * 



(4) 



where 



\hL\' 

t — m?. 



(5) 



and M is the invariant mass of the l^l pair. The QCD corrections to the Drell-Yan 



process have been calculated [|17| to the next-to-leading order and is a function of the 
cm. energy (i/i) of the collider, the structure functions used and the subprocess scale 
M . For M^20 GeV, the dependence on M is marginal and one may approximate it by 
a scale-independent constant. In the absence of an explicit calculation of the higher- 
order effects, we assume that the QCD correction for the t-channel scalar exchange 
process is the same as that for the Drell-Yan process. One may argue that in the 
presence of scalars, these corrections are likely to be even larger. 

To obtain the cross section for a hadronic collision AB l^l^X, we have to 
convolute the expression in eq. (H) with the parton distributions. For example, the 
invariant mass distribution is given by 



da 2M 



dM 



J dxdt 



qA{x)qB{Mysx)^ + {A 
dt 



B] 



(6) 



where x is the fraction of the momentum of A carried by q^. While the sum in eq.(^ 
runs over all species of quarks for the SM piece, for the scalar contribution it runs only 
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over those quarks with which the scalar couples. For the rest of this analysis, we shall 
concentrate only on those that couple to one or both of the first generation quarks. 



We use the MRSD-' structure functions |18] evaluated at the scale Q = M . For this 



choice, the QCD isT-factor ~ 1.3(1.1) for ^ = 1.8(14) TeV. 

The CDF collaboration at Fermilab measures the Drell-Yan cross section [IB 
da/dM in the range 11 < M < 150 GeV. They restrict their measurements to |?7±| < 1, 
where r/± are the pseudorapidities of respectively. Furthermore, the numbers quoted 
are averaged over the rapidity of the lepton pair. For ease of comparison, we shall 
consider the same distribution (including the actual binning for M as used in the ex- 
periment). In Fig. we compare the experimental data with the theoretical curves 
for the SM, and those for scalars of Type IVa (see Table We have set here = 1, 
i.e. the Yukawa coupling is of the electroweak strength. Though the effect of the 
scalar starts becoming visible at around M ~ 65 GeV, it is immediately swamped by 
the Z-pole. At large M values though, the deviation is significant and information 
about scalars may be extracted. However, as even a cursory glance at the figure would 
reveal, the data at present are rather poor and any quantitative statement would be 
premature. We rather wish to point out that a refinement of measurement at the 
high invariant mass end (M ^ 150 GeV) of the spectrum (along with an increase of 
luminosity at the possible Tevatron upgrade) would enhance considerably the ability 
to detect such scalar particles. 

We now turn to the LHC. Due to the large operating energy (14 TeV) one would 
expect to see a more pronounced effect. In our analysis, we assume an integrated 
luminosity of 100 fb~^. In Fig. 0, we first compare the invariant mass distribution for 
different masses of a Type IVa scalar [h^ = 1) with that for the SM. We use kinematic 
cuts of \rj±\ < 3, numbers that have often been quoted. As is expected, the effect of 
the scalar exchange is more pronounced for large M. Due to the much larger energy 
available (and partly due to the wider angular coverage), the presence of much heavier 
scalars can be detected. 

The information in Fig. ^ can be used to put bounds on the two-parameter space 
(m^, h) — here h = h^^R, as the case may be — that can be achieved at the LHC. 
It is obvious that imposing an acceptability cut on the invariant mass would enhance 
the effect of the scalar. In our analysis, we arbitrarily set Mmin = 500 GeV. Denoting 
the number of events expected in the SM, and in a theory with such a scalar by usm 
and respectively, we demand that 

1^0 - nsM\ > ^V^SM (7) 

for the effect to be considered visible. This allows us to draw contours in the m^-h 

■^It should be noted that the sharpness of the Z-resonance is somewhat reduced by the rather 
coarse binning adopted in the experiment. 
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plane (Fig. ^ for various types of scalars. The region of the parameter space above 
the respective curves can then be ruled out at the 3cr level. The difference in the 
sensitivity to the scalar type is a consequence of the difference in their coupling, as 
also of the relative abundance of various quarks in the proton. While it is true that a 
2.5 TeV scalar may easily be pair produced at the LHC, such an event would afford 
us no handle on its couplings to the SM fermions. Indeed, the bounds presented here 
are the strongest that can be achieved on leptoquark couplings in the near future. 
For the I/ip case, this is illustrated even better. As we have pointed out right at the 
beginning, the low energy constraints on ]/ip couplings are stronger than those for 
the usual leptoquarks. For our interaction, is phenomenologically identical to a 
case where both Type II a and Type III scalars are present, with the contributions 
adding incoherently. If we assume the two squarks {ul and dji) to be mass-degenerate, 
the bound on the coupling lies in between the two above curves, and is somewhat 
weaker than those derived from low energy processes One should realise though 
that this experiment provides an independent constraint on the parameter space for the 
theory. Interestingly, in the extreme case of m{dfi) ^ m{uL), the low energy constraint 
becomes relatively unimportant and the experiment discussed here would provide the 
strongest constraints. Indeed, for all such scalars except that of Type ///, dilepton 
data at the LHC would lead to bounds significantly stronger than those obtained from 
low energy processes such as vr-decay, charge-current-universality etc. 

As we had mentioned, our choice of Mmin was rather arbitrary. The exact value 
that should be adopted to maximize the effect is dependent on the value of that one 
aims to probe. On the face of it, a larger value of Mmin should serve to eliminate more 
of the SM contribution. This advantage can, however, be nullified by the consequent 
loss of statistics. In fact, a larger Mmin is more useful for exploring larger m,^. To 
illustrate our point, we superimpose the 3a contours for two different values of Mmin 
in Fig. 0. 

Until now, we have neglected the other interesting dynamical variable in the pro- 
cess, namely the difference [At] = rjj^—rjJ) of the lepton rapidities. Since this is nothing 
but the scattering angle in the subprocess cm. frame, one expects to see the difference 
between an s-channel and a t-channel process. This is borne out strikingly in Fig. ^. 
We would like to point out that instead of comparing the integrated cross section, 
much more information can be gleaned from a study of the differential distribution 
(Pa/dM dArj. This can be done, at a quantitative level, by performing a test. We 
refrain from doing so, however, as such an analysis presumes some knowledge of the 
detector. We prefer instead, to only point out that once the detector parameters are 
well understood, such an analysis should be undertaken so as to improve the sensitivity 
of the experiment. 

In summary, we have examined the virtual effect of a t-channel scalar exchange 
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in dielectron production at CDF and tried to foresee its possible impact at the LHC. 
In contrast to the pair production mechanism, the virtual effects are very sensitive to 
its Yukawa coupling. We point out that a few precise measurements of the dilepton 
distribution at CDF at higher invariant masses (M ~ 150 GeV) would allow us to 
probe the existence of a scalar of mass upto a few hundred GeV and couphng the 
electron and a first generation quark with electroweak strength. Similarly, at the 
LHC, we observe the possibility of probing scalar masses upto (1-3) TcV for a Yukawa 
coupling of electroweak strength. These estimates are significantly better than the 
bounds inferred from analyses of low-energy processes. We would also like to point 
out that such an analysis can as easily be done for a dimuon pair, and the bounds 
would be similar. For a t'^t~ pair, the situation would be a bit more complicated as 
T-identification efficiency is likely to be lower. If this efficiency could be raised, such 
an experiment would lead to a significant improvement in constraints on such scalars 
coupling a r to a quark. Finally, one may also investigate possible couplings of leptons 
to the heavier quarks, although these bounds would be weaker on account of the lower 
densities within the proton. 
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(2, -7/6) 
(1 1/3) 




II a) 
h) 


eR^dRcf) 


(2, -1/6) 
(1, 4/3) 


Corresponds to ul 


III 




(1, 1/3) 


Corresponds to dR 


IV a) 
b) 




(2, -7/6) 
(3, 1/3) 




V a) 
b) 


GRQlft) 

Il''Ql4' 


(2, -7/6) 
(3, 1/3) 


m(0^/^) > m(0^/^) 



Table 1: The possible scalar couplings between a charged lepton and a quark, grouped 
according to their chiral structures. For Type IV, we assume all the scalars to be mass 
degenerate. Type V corresponds to the (unlikely !) case where only the coupling to 
di is of importance. 
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Figure 1: The dilepton invariant mass distribution at CDF. The curves correspond 
to the theoretical expectation for the SM and for leptoquarks (Type IVa) of masses 
= 100, 200, 300 GeV respectively. For M < 150 GeV, the invariant mass binning 
is the same as that used in . 
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Figure 2: The dilepton invariant mass distribution at the LHC The curves correspond 
to the theoretical expectation for the SM and for Type IVa scalars of masses m<^ = 
500, 1500, 2500 GeV respectively. 
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Figure 3: The parameter space that can be probed at the LHC The different curves 
are for the various scalar types hsted in Table |l], with the ordinate corresponding to 
the respective coupling. The part of the parameter space above the individual curves 
can be ruled out at the 3a level. 
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Figure 4: The dependence of the exploring abihty on the invariant mass cut. The two 
curves (Sa) shown are for the Type IVa scalar (see Table P and for Mmin = 500 GeV 
and 1 TeV respectively. 
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Figure 5: The differential cross section for various values of scalar masses (Type IVa) 
as a function of the difference of the lepton rapidities. 
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